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Abstract 

We describe the linear and nonlinear stabilities of certain configurations of point vortices on 
the sphere forming relative equilibria. These configurations consist of up to two rings, with and 
without polar vortices. Such configurations have dihedral symmetry, and the symmetry is used to 
block diagonalize the relevant matrices, enabling their eigenvalues to be calculated. 

1 Introduction 

Since the work of Helmholtz [ H | systems of point vortices on the plane have been studied as finite- 
dimensional approximations to vorticity dynamics in ideal fluids. For a general survey of patterns of 
point vortices see [ANSTV|. Point vortex systems on the sphere, introduced by Bogomolov [B77|, 
provide simple models for the dynamics of concentrated regions of vorticity, such as cyclones and hur- 
ricanes, in planetary atmospheres. In this paper we consider a non-rotating sphere, since the rotation 
of the sphere induces a non-uniform background vorticity which makes the whole system infinite- 
dimensional. 

As in the planar case, the equations governing the motion of N point vortices on a sphere are 
Hamiltonian IB77I and this property has been used to study them from a number of different view- 
points. Phase space reduction shows that the three vortex problem is completely integrable on both 
the plane and the sphere: the motion of three vortices of arbitrary vorticity on a sphere is studied in 
[KN98|. The stabilities of some of the relative equilibria given in [KN98| are computed in [PM98| 
and numerical simulations are presented in [MPS99]. The existence of relative equilibria of A" vor- 
tices is treated in [LMR01|, and the nonlinear stability of a latitudinal ring of A" identical vorticities 
is computed in [BC03|, and independently in the present paper. It follows that the linear stability 
results obtained by [ PD93 1 are in fact Lyapunov stability results. The stability of a ring of vortices on 
the sphere together with a central polar vortex is studied in [CMS03|, and again independently in the 
present paper. The existence and nonlinear stability of relative equilibria of A" vortices of vorticity +1 
together with A" vortices of vorticity —1 are studied in [LP02|. It has also been proved in [LP| that 
relative equilibria formed of latitudinal rings of identical vortices for the non-rotating sphere persist 
to relative equilibria when the sphere rotates. Finally, studies of periodic orbits of point vortices on 
the sphere can be found in ISTIIToOlllLPthl . 

Our study of the stability of relative equilibria is based on the symmetries of the system, and 
especially the isotropy subgroups of the relative equilibria. The Hamiltonian is invariant under ro- 
tations and reflections of the sphere and permutations of identical vortices. However, some of these 
symmetries (eg reflections) are not symmetries of the equations of motion: they are time-reversing 
symmetries. From Noether's theorem, the rotational symmetry provides three conserved quantities, 
the components of the momentum map $ : V — > R 3 where V is the phase space. 
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Relative equilibria are dynamical trajectories that are generated by the action of a 1 -parameter 
subgroup of the symmetry group. More intuitively, they correspond here to motions of the point 
vortices which are stationary in a steadily rotating frame. In other words, the motion of a relative 
equilibrium corresponds to a rigid rotation of N point vortices about some axis (which we always 
take to be the z-axis). In the same way as equilibria are critical points of the Hamiltonian H, relative 
equilibria are critical points of the restrictions of H to the level sets <J> -1 (//). Section|2]is devoted to 
a description of the system of point vortices on the sphere, and to an outline of stability theory for 
relative equilibria. The appropriate concept of stability for relative equilibria of Hamiltonian system 
is Lyapunov stability modulo a subgroup. The stability study is realized using on one hand the energy 
momentum method [Pa92 Or98 ] which consists of computing the eigenvalues of a certain Hessian, 
and leads to nonlinear stability results, and on the other hand a linear study computing the eigenvalues 
of the linearization of the equations of motion. To this end, we block diagonalize these matrices using 
a suitable basis, the symmetry adapted basis (Section which makes use of the specific dihedral 
isotropy subgroup of the relative equilibrium. 

The remaining five sections each treat one of five different types of relative equilibria, consisting of 
rings of identical vortices together with possible vortices at the poles, whose existence were proved in 
[LMR01 1. The notation for the different configurations is taken from the same source and is described 
at the end of the introduction. We now outline the main stability results. 

We begin in Section |4] by computing the stability of the relative equilibria consisting of a single 
ring of identical vortices, a configuration denoted C nv (R) (Figure [TTTT a)). We show in Theorem 14.21 
that for n > 7, they are unstable for all co-latitudes of the ring, while for n < 7 there exist ranges of 
Lyapunov stability when the ring is near a pole. These results are not new IPD931 IBC03I . but serve to 
demonstrate the method used in later sections. 

In Section |5j we study the stability of the relative equilibria C nv (R,p) (Figure [TTTT b)) which are 
configurations formed of a ring of n identical vortices together with a polar vortex. For n > 7 they 
are all unstable if the vorticity k of the polar vortex has opposite sign to that of the ring. However 
if the vorticities have the same sign then for each co-latitude of the ring there exists a range of k for 
which the relative equilibrium is Lyapunov stable. Adding polar vortices can therefore stabilize the 
unstable pure ring relative equilibria. The detailed results are contained in Theorem l5.2l its corollary 
and the following discussion. Our results are consistent with those of [CMS03|, though our methods 
give more regions of stability than theirs — see Remark \5A\ 

In Section |6] we obtain analytic (in)stability criteria for the relative equilibria C nv (R, 2p) which 
are configurations formed of a ring of n identical vortices together with two polar vortices, but only 
with respect to certain modes (£ > 2). As in the case of a single polar vortex, the two polar vortices 
play the role of control parameters for the stability. The details are contained in Theorem 16.21 A 
numerical investigation is needed for the remaining (£ = 1) mode in order to provide stability criteria; 
this is being pursued separately. 

Finally, in Sections and [8] we investigate configurations formed of two rings of arbitrary vor- 
ticities (each ring, as always, consisting of identical vortices). In [LMR01 1 it was shown that two 
rings of n vortices can be relative equilibria if and only if they are either aligned or staggered. These 
two arrangements are denoted C nv (2R) and C nv (R, R r ) respectively (see Figure [L2l . Here we show 
that for almost all pairs of ring latitudes there is a unique ratio of the ring vorticities for which these 
configurations are relative equilibria. Numerical computations of their stabilities suggest that these 
relative equilibria can only be stable if n < 6, and in the aligned case the two rings must be close to 
opposite poles, and hence have opposite vorticities. In some cases, staggered rings may also be stable 
when in the same hemisphere. 
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(a)C iv (2R) (b)C 4 v(R,R') 
Figure 1.2: Configurations of types C nv (2R) (2 aligned rings) and C nv (R, R r ) (2 staggered rings). 
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Symmetry group notation All possible symmetry types of configurations of point vortices on the 
sphere were classified in [LMR01 1. The symmetry group of the system is of the form 0(3) x S, where 
S is a group of permutations, and a particular configuration with symmetry, or isotropy, subgroup 
£ < 0(3) x S is denoted T(A), where T is the projection of S to 0(3) and A represents the way S 
permutes the point vortices. The classical Schdnflies-Eyring notation for subgroups of 0(3) is used. 

In this paper we single out configurations consisting of concentric rings of identical vortices, with 
the same number of vortices in each ring, and with possible polar vortices. These configurations have 
cyclic symmetry (in the "horizontal plane"), and the Schonflies-Eyring notation for this subgroup of 
0(3) is C n < SO (3). In fact we only consider the cases where the rings are either aligned (the 
vortices lie on the same longitudes) or staggered (they lie on intermediate longitudes, out of phase 
by 7r/n). In this case the symmetry group is the larger dihedral group C nv (n being the number of 
vortices in each ring, and v denoting the fact that there are vertical planes of reflection). For such 
configurations, we write C nv {k\R, &2-R', k p p) to mean that there are k r = k\ + &2 rings and k p polar 
vortices. The difference between R and R' is that the k\ rings R are aligned and the k2 rings R' are 
staggered with respect to the first (and so aligned with each other). Of course k p = 0, 1 or 2. 

2 Point vortices on the sphere and stability theory 

In this section we briefly recall that the system of point vortices on a sphere is a n-body Hamiltonian 
system with symmetry and we review the stability theory for relative equilibria. 

2.1 Point vortices 

Consider n point vortices xi,...,x n € S 2 with vorticities k%, . . . , n n € R. Let 0j, fa be respec- 
tively the co-latitude and the longitude of the vortex X{. The dynamical system is Hamiltonian with 
Hamiltonian given by 

H = — K i K j m (l ~~ cos Qi cos &j ~ sm Qi srn 8j cos(0j — 4>j)) 

and conjugate variables given by qi = y/\Ki \ cos 0j and Pi = sign(Kj)y / |Kj|(/)j. 

The phase space is V = {(x±, . . . , x n ) € S 2 x • • • x S 2 \ xi ^ Xj if i ^ j} endowed 
with the symplectic form uj = ^ Kjsin6*j dOi A dfa. The Hamiltonian vector field Xh satisfies 
uj( ■ , Xh(x)) = dH x . If we consider S 2 as a subset of ]R 3 , so the vortices xj G M 3 , then we obtain 

x,i = X H (x)i = } kj — , i = l,...,N, 

H = — KjKj ln(||xj — Xj\\ 2 /2). 

i<j 

It follows that H is invariant under the action of 0(3). The symplectic form is S0(3) invariant and 
so Xh is S0(3)-equivariant. The reflections in 0(3) reverse the sign of the symplectic form and so 
are time-reversing symmetries of Xh- Moreover H, to and Xh are all invariant or equivariant with 
respect to permutations of vortices with equal vorticity. 

The rotational symmetry implies the existence of a momentum map $:'P^so(3)*~]R 3 : 

N 
3=1 
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which is conserved under the dynamics. In other words each of the three components of is a 
conserved quantity. 

2.2 Relative equilibria 

A point x e G V is a relative equilibrium if and only if there exists £ G so (3) ~ R 3 (the angular 
velocity) such that x e is a critical point of the function Hg(x) = H(x) — (&(x), £), where the pairing 
( , ) between R 3 and its dual is identified with the canonical scalar product on R 3 . Equivalently, 
relative equilibria are critical points of the restriction of H to $ _1 (//), since the level set Q~ l (p) are 
always non-singular for point vortex systems of more than two vortices. The function is called the 
augmented Hamiltonian. 

Since the momentum is conserved, we can choose a frame for R 3 such that is parallel to the 
z-axis (provided the momentum is non-zero). It follows from the symmetry that the angular velocity 
£ € R 3 is also parallel to the z-axis. We can therefore identify £ and <3? with their z-components and 
the augmented Hamiltonian becomes simply H^(x) = H(x) — £$(x). 

Let / : V — ► R be a iv~-invariant function with K a compact group. Recall that Fix(iv") = {x G 
V | g ■ x = x, \/g G K}. The Principle of Symmetric Criticality jP79| states that a critical point 
of the restriction of a iT-invariant function / to Fix(K) is a critical point of /. As a corollary, if the 
Hamiltonian is invariant under K and x e is an isolated point in Fix(ET) n <&~ l {p), then x e is a relative 
equilibrium. It follows in particular that all configurations of type C nv (R), C nv (R,p) (Figure fT7TT >. 
and C nv (R, 2p) are relative equilibria: take K such that tt(K) = C nv where it : 0(3) x S n — ► 0(3) 
is the cartesian projection. 

Finally, one can show that if x e is a relative equilibrium with angular velocity £, then H% is a 
G Xe -invariant function. 

2.3 Stability theory 

Stability is determined by the energy momentum method together with an isotypic decomposition of 
the symplectic slice. We recall the main points of the method. 

Let x e G V be a relative equilibrium, p, = $(x £ ), and £ its angular velocity. The energy momen- 
tum method consists of determining the symplectic slice 

N = (so(3) M • Xe) 1 - DKerD^(x e ) 

transversal to so(3)^ • x e , where 

S0(3)^ = {g G S0(3) | Coad 9 -p = p} 

and then examining the definiteness of the restriction d 2 H^\j^(x e ) of the Hessian d 2 H^(x e ) to J\f. If 
K is a group acting on the phase space a relative equilibrium x e is said to be Lyapunov stable modulo 
K if for all if-invariant open neighbourhoods V of K • x e there is an open neighbourhood U C V 
of x e which is invariant under the Hamiltonian dynamics. The energy momentum theorem of Patrick 
IPa92 l holds since S0(3) is compact, and so we have: 

If d 2 H^\j\f{x e ) is definite, then x e is Lyapunov stable modulo S0(3)^. 

For p ^= 0, SO(3) M is the set of rotations with axis (p), and so isomorphic to SO(2), while for 
p = 0, SO(3) /t = SO(3). If p Lyapunov stability modulo SO(3) M of a relative equilibrium with 
non-zero angular velocity corresponds to the ordinary stability of the corresponding periodic orbit. 
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The second step consists of performing an isotypic decomposition of the symplectic slice J\f in 
order to block diagonalize d 2 H^\j^(x e ). Let V be a finite dimensional representation of a compact 
Lie group K. Recall that a K invariant subspace W C V of K is said to be irreducible if W 
has no proper K invariant subspaces. Since K is compact, V can be expressed as a direct sum of 
irreducible representations: V = W\ © • • • © W n . In general this is not unique. There are a finite 
number of isomorphism classes of irreducible representations of K in V, say U\, . . . , U%. Let Vf. 
(k = 1, ...,£) be the sum of all irreducible representations Wj C V such that Wj is isomorphic 
to Uf.. Then V = V\ © • • • © Vj>. This decomposition of V is unique and is called the K -isotypic 
decomposition of V [Se78 ]. By Schur's Lemma, the matrix of a K-equivariant linear map / : V — > V 
block diagonalizes with respect to a basis B = {B\, . . . , B{\ where Bk is a basis of Vfc, each block 
corresponding to a subspace Vfe. The basis £> is called a symmetry adapted basis. 

Let G denote the group of all symmetries of H and Xh and G x the subgroup consisting of time- 
preserving symmetries. In the case of N identical vortices we have G = 0(3) x SV and G x = 
S0(3) x Sn- Since is a G Xe -invariant function, d 2 H^\j^(x e ) is G Xe -equivariant as a matrix. 
Moreover the symplectic slice J\f is a G Xe -invariant subspace and so we can implement a G Xe -isotypic 
decomposition of M to block diagonalize d 2 H^\j^(x e ). This block diagonalization of d?H^\j^{x e ) 
simplifies the computation of its eigenvalues, and hence its definiteness. If it is definite then the 
relative equilibrium is Lyapunov stability modulo S0(3)^. 

If d 2 H^\j^f(x e ) is not definite then we study the spectral stability of x e . In particular we ex- 
amine the eigenvalues of Lj^, the matrix of the linearized system on the symplectic slice, that is 
Lj\f = 3j^d 2 H^\_\f(x e ), where J^ 1 is the matrix of w|jv- The matrix Ljv" is G Xe -equivariant and so 
we perform a G Xe -isotypic decomposition of M to obtain a block diagonalization of Lj^f, and so to 
determine the spectral stability of x e . In particular, if Lj\f has eigenvalues with non-zero real part, 
then x e is linearly unstable. Note that the block diagonalization of d 2 H^\j^(x e ) refines that of Lj^f 
since G£ e C G Xe . 

Throughout this paper, Lyapunov stable will mean Lyapunov stable modulo S0(2) unless specified 
otherwise, and p = $(x e ) will denote the momentum of the configuration x e . We will assume in most 
of the results that the relative equilibrium has a non-zero momentum. It is straightforward to check that 
almost all points of type C nv (R, k p p) have a non-zero momentum, so this is not a strong assumption 
for that case. This statement is also true for C n „(2i?) and C nv (R, R') relative equilibria, as we shall 
see. 

3 Symmetry adapted bases for rings and poles 

In this section we give the ingredients needed to determine the symmetry adapted bases for the sym- 
plectic slice at the configurations described above, that is those of type C nv (k\R, k%R! ', k p p). In the 
first subsection we give a general symmetry adapted basis for the tangent space T Xe V to the phase 
space at such a configuration, and express the derivative of the momentum map and tangent space to 
the group orbit in this basis. In the following two subsections we describe the isotypic decomposi- 
tion of T Xe V, first for a single ring and then in general. Recall that the isotropy subgroup is always 
a dihedral group C nv and that the irreducible representations of this group are of dimension 1 or 2. 
The actual symmetry adapted bases of the symplectic slices will be given case-by-case in the follow- 
ing sections. We do not give the proof of the results in the first subsection, since they can be easily 
deduced from the proofs of Propositions 4.1, 4.2, 4.3, 4.4 in ILP02I . 
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3.1 Description of the symplectic slice 

Let x e be a C nv (k\R, k 2 Rf, k p p) configuration. The total number of vortices is N = n[k\ + k 2 ) + k p . 
We suppose the vorticities in the k r = k\ + k 2 rings are Kj for j = 1, . . . , k r while the vorticities of 
the possible polar vortices are At at for the North pole and k$ for the South pole. 

For each ring j = 1, . . . , k r let s = 1, . . . , n label the vortices in the ring in cyclic order, and 
define tangent vectors in T Xe V by 



where £ = 0, . . . ,n — l,i = and (jp, = or ir/n depending on whether the j th ring of vortices is 



3 

of type R or R'. Note that j3fg , (3^ vanish for i = and n/2 (for n even). For each pole j = 1, . . . , k. 
we also have tangent vectors 6xj and Syj . 

The tangent vectors defined in the last paragraph are almost canonical, in the sense that 



if I = 0,n/2 



W3#) {'-^P-Kj otherwise 

= 

u(Sxj,6yj) = sign(zj)Kj, 

while the other pairings vanish. 

In order to compute a specific basis for the symplectic slice it is necessary to have expressions for 
the derivative of the momentum map and the tangent space to the group orbit at a C nv (k\R, k 2 R' ', k p p) 
configuration. These expressions are given in the next two propositions. Since C nv refers to fixed 
vertical reflection planes, the values of 0° in (I3.lt above can be taken to be: 



if i = 1 . . . fci 

tt/u if j = (k% + 1) . . . k r 



where k r = k\ + k 2 is the total number of rings. 

Proposition 3.1 At a C nv (k\R, k 2 R', k p p) configuration, the differential of the momentum map is 
given by, 

d$>(5x, Sy, a, /?) 

= Y K i ( 5x i + i6y i ) + Y K i cos 6 i ( a S + ) + * Y K i sin ^' ( a S + ) 
j polar i =1 i =1 



■ (o) 



where the direct sum corresponds to the C nv -invariant decomposition of so (3)* as a direct sum of a 
plane and the line Fix(C n t,,so(3)*) (the "z-axis"). 
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Proposition 3.2 Let x e be a C nv (kiR, &2-R', k p p) configuration, and p = $(x e ). If p ^ 0, then the 
tangent space to the orbit so(3) M • x e is generated by the vector 

If p = 0, then so (3)^ • x e = so (3) • x e and is generated by the three vectors: 

£4 



i- 1 jpolar 



w/iere o«e mM5t fa&e ?/ze rea/ araci imaginary parts of the first line. 
3.2 A single ring 

Let x e be a C nv (R, k p p) configuration, that is a single ring together with k p polar vortices where, 
of course, k p = 0,1 or 2. Since there is only one ring, we write offi , ffi , , instead of 
af^g, f3fl,a>f^, (3^ in order to lighten the formulae. 

The irreducible representations (subspaces) of dimension 1 of the action of G Xe on T Xr V are: 



«i° / j \ a i / if n is odd, 



while they are 



(0)\ / (0)\ / (n/2)\ / (n/2)\ ., . 

a 9 /> \ a 4> /' \ a e /> \ a V / if n is even. 



Of these, only (a^ 2 ^ and (a^ 2 ^ He in the symplectic slice. Moreover, (cxg 1 ^ 2 ^ and (ot^ 2 ^ 1 
are G Xe -isomorphic representations, but not G Xe -isomorphic representations. 



The following spaces are irreducible representations of dimension 2 of the action of G Xe on T Xe V: 
[af , f3f ] ) , (af ,pf),l<t<n-\,i£ n/2, {Sxj,6 yj ) ,j = l,...,k p . 

The representations (o^\ > ( a ^'/^^) > i^j^Vj) do not lie in the symplectic slice, while 

the others do. Moreover, (otg\ fijp^, (a^\f3^^, and (5xj,5yj) are G Xe -isomorphic representa- 
tions. 

3.3 General case 

In the general case where x e is a C nv (k\R, k^R! ', k p p) configuration (k p = 0, 1 or 2), we have the 
following decomposition. 
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The subspaces {or- j j , ^aA ^ \ , j = 1, . . . , k r are G Xe -irreducible representations of dimension 1 

and are Gx e -isomorphic representations. If n is even, then we have 2k r additional G Xe -irreducible 
representations 

m<i =*<■■■<'<' 

which are Gx e -isomorphic. 

The subspaces (uffiiPffi) > ( a f]f» > j = ^ - • • > &r, I < £ < n - l,£ ^ n/2 and (Sx r , Sy r ) , 
r = 1, . . . , k p are G Xe -it reducible representations of dimension 2. Moreover, the subspaces 

(<*?}>$)> (4>S)> ( s *r,sy r ), 

(j = l,...,k r , r = 1, . . . , k p ) are Gx e -isomorphic representations. 

The difficulty in the case of several rings is that usually the subspaces listed above do not lie in the 
symplectic slice. One needs therefore to find linear combinations of the above vectors that do lie in 
the symplectic slice, and that the irreducible and isomorphism properties of the representations are 
preserved. We will give the symmetry adapted bases for the cases of two rings without polar vortices, 
that is C nv (2R) and C nv (R, R') configurations, in the last two sections. 



4 A ring of identical vortices: C nv (R) 

The linear stability of C n v(R) relative equilibria was determined by Polvani and Dritschel in I PD93I . 
Recently, Boatto and Cabral [BC03] studied their Lyapunov stability and found that the two types of 
stability coincide: whenever the relative equilibrium fails to be Lyapunov stable the linearization of 
Xh has real eigenvalues. In this section, we give another proof using the geometric method of this 
paper. 

For n vortices of unit vorticity the Hamiltonian is 

H(9j,(pj) = — In (1 — sin 9j sin 6k cos(cpj — cfrk) — cos 9j cos 6k) 

j<k 

and the augmented Hamiltonian is = H — £ ]TV cos 6j. 

Let x e be a C nv (R) relative equilibrium and 6 the co-latitude of the ring. The angular velocity of 

x e is 

(n — 1) cos #o 
* = sffi^ ' 

since has a critical point there and ^^{x e ) = ( n ~ 1 ) c ° s ^~? sm 
The second derivatives of H at the relative equilibrium are: 

d 2 H (n-l)(n-5) d 2 H 1 

86'j 6sin 2 6)o d0jde k 2 sin 2 6 sin 2 (7r(j-fc)/n) 

d 2 H _ o d 2 H _ n 

d 2 H _ 1 d 2 H _ 1 

d(j>'j Z^r=l 2sin 2 (r7r/n) dc/>jd(p k 2 sin 2 -k)/n) ' 

We note that Y,%1 l/sm 2 (vrr/n) = \(n 2 - 1) IH751 . 
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Notation. In order to harmonize the statements of the results between n even and n odd, we intro- 
duce the following notation: let T][ , in^ > V3 j V4 t> e objects defined for all 2 < I < [^^p] , where 
[m] is the integer part of m G N, and only 7/^ , 773 for ^ = n/2 when n is even. Then define 

{r/f ) ,r?f ) I 2 < £ < [n/2]}* to be 



or even n : {n? >, rg> , r,?\ | 2 < * < § - l} U ^} (4 ^ 

for odd n : {r?f } , 7$ } , ^ , r/f } | 2 < * < [n/2] } . 



The following proposition gives the symmetry adapted basis for a C nv (R) configuration. 
Proposition 4.1 Assume ^ / 0. With respect to the following basis for the symplectic slice M at x e , 

( ei ,e 2 , {aWaW^,^)|2</<[n/2]}*) 

where 

e\ = sin 9q o>q + cos 9q 
e 2 = sin 9o pp — cos 6$ 

f/je Hessian d 2 H^\j^(x e ) block diagonalizes in 1 x 1 blocks, and Ljsf block diagonalizes in 2 x 2 

Proof. It is straightforward to check that the vectors above do form a basis for the symplectic slice 
at x e thanks to Propositions 13. Il and l3.2l 

The Hessian d 2 H^\j^{x e ) and the linearization Lj^f are both G% e -invariant. Assume n odd. It fol- 
lows from Section l3~2l and Schur's Lemma (see the Introduction) that d 2 H^\j\f(x e ) and L/j- both 
block diagonalize into 4x4 blocks and one 2x2 block corresponding to the subspaces Vi = 
(affif3^\a^,(3^j and (ex, e 2 ), respectively. See the proof of Theorem 4.5 of ILP02I for a de- 
tailed proof of a similar assertion. 

Now fix I and denote by s an anti-symplectic (time-reversing) element of G Xe . For example s 
could be the reflection y ^ —y together with an order two permutation of S n . The restriction of 

to Vi is Z 2 [s] -invariant. Moreover (<4 >/^i / anc ^ \P(P i a 6/ WQ n on-isomorphic irreducible rep- 



V / 

resentation of Z 2 [s] on Vg. Hence d 2 H^\j s f(x e ) block diagonalizes into 2x2 blocks which correspond 
to subspaces (o£q \ an d (ei, e 2 ). This result does not depend on the details of the 

Hamiltonian, only its symmetries. However taking in account its particular form, one can improve the 
block diagonalization. Indeed one has 

d 2 Hi{x e )-{af,pf) = d 2 H^x e )-{f3f\af) = 

and d 2 H^(x e ) • (ei, e 2 ) = which gives the desired diagonalization of the Hessian. 

The particular form of the symplectic form also enables us to improve the diagonalization of Lj^. 
Among the basis vectors of Vi, only u(ay , affl ) and , fijp ) do not vanish, and so the restriction 

/ (£) (£) 

of u) to Ve block diagonalizes into two 2x2 blocks which correspond to the subspaces (a\ ,a\ 



and The block diagonalization of Lj\f then follows from Ljv = J^rd? Hg\j\f(x e ). 

The case n even is very similar, except that there is an additional 2x2 block in the Gx e -isotypic 
decomposition, and leads to the same result. □ 
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The block diagonalization of d 2 H^\j\f(x e ) and La/ enable us to find formulae for their eigenvalues, 
and thus to conclude criteria for both Lyapunov and linear stability. 

Theorem 4.2 The stability of a ring of n identical vortices depends on n and the co-latitude 9q as 
follows: 

n=2 is Lyapunov stable at all latitudes; 
n=3 is Lyapunov stable at all latitudes; 

n=4 is Lyapunov stable if cos 2 Oq > 1/3, and linearly unstable if the inequality is reversed; 
n=5 is Lyapunov stable if cos 2 9q > 1/2, and linearly unstable if the inequality is reversed; 
n=6 is Lyapunov stable if cos 2 0q > 4/5, and linearly unstable if the inequality is reversed; 
n>7 is always (linearly) unstable. 

Proof. Any arrangement of two vortices is a relative equilibrium [KN98 1. When perturbing such a 
relative equilibrium, we obtain a new relative equilibrium close to the first. Thus any relative equilib- 
rium of two vortices is Lyapunov stable modulo SO(2) (modulo SO(3) if p = 0). 

Hence assume n > 3. We first study Lyapunov stability. Suppose further that / and so the 
ring is not equatorial. A simple calculation shows that d 2 H^(x e ) • , pp) = d 2 H^(x e )-(af\af^) 
and d 2 H^(x e ) ■ {/3^\p { ^) = d 2 H^(x e ) ■ {af ,af). Hence it follows from Proposition O that 

d 2 HzW(x e ) = diag (a^Ax.IA^A^.A^A^ | 2 < I < [n/2]}*) 
(recall notation from (14.11) ) where 

Ai = sin^oA^ +cos 2 6 X ( ^, 
\ e — a H^{x e ) ■ (a e ,a e J, 
\<f> = d 2 H,(x e ).(af,af). 
Thanks to the following formula [H73J (p-271) 

sm (Trj/n) o 
we find after some computations that A^ = n£(n — tj/2 and 

X f = H £ - l)(n - £ - 1) + (n - 1) cos 2 flo] . 

z sm uq 

The eigenvalues A^ are all positive and Ai = n(n — 1) cos 2 6q > 0, thus the relative equilibrium is 
Lyapunov stable (modulo SO(2)) if (n - 1) cos 2 9 > (£- l)(n - 1- 1) for all I = 2, . . . , [n/2], that 
is if cos 2 #o > ([n/2] — l)(n — [n/2] — l)/(n — 1). This gives the desired values. 
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We now turn to linear stability. It follows from Proposition 14.11 and the block diagonalization of 
d 2 H^(x e ) that 




-A?> \ / -A<« 



where the blocks are given up to a strictly positive scalar factor. The eigenvalues of Lj^ are therefore 



±i\xA±iyJ\f\f \ 2<£<[n/2}^ 

(£) 

(up to a positive factor) and so the relative equilibrium is linearly unstable if X e > for some £, that 
is if 



cos 2 9q < 



n- 1 V L 2 



n 



1 



n + 1 



1 



In particular this inequality is satisfied if 6q = it/ 2 and n > 3. 

When the ring is equatorial, one has 8q = tt/2 and fi = 0. In particular Ai = 0. This is because the 
symplectic slice is smaller (G^o = SO(3)): it follows from Proposition 13.21 that we have to remove 
the vectors e\, e2 from the basis for fi ^ (that is to remove Ai from the previous eigenvalue study). 
It follows that the C nv equatorial relative equilibria are linearly unstable for n > 3, and Lyapunov 
stable (modulo SO(3)) for n = 3. □ 

The proof shows that the 'critical mode' for stability is £ = [n/2]. For n > 7 a ring is always 
unstable to this mode, while for 4 < n < 6 the ring first loses stability to this mode as it moves closer 
to the equator. This loss of stability is accompanied by a pitchfork bifurcation to a pair of staggered 
rings when n = 4 or 6 (i.e. to types C2 V (R, R') and C% V (R,R') respectively, in the notation of 
[LMR01 1). In the case n = 5 the bifurcation is transcritical to an "equatorial" vortex and two pairs 
that are reflections of each other in that equator (ie type C/ l (2i?, E) in [LMR01 1). The bifurcations 
for n = 4 and 5 are illustrated in Figures 7 and 8 of [LMR01 1 respectively. 

5 A ring and a polar vortex: C nv (R, p) 

We assume that the polar vortex lies at the North pole and its vorticity is k, while the other vortices 
are all identical with vortex strength 1. The relative equilibrium is of symmetry type C nv (R,p) and 
denoted x e . In this case, the Hamiltonian is given by 

H = H r + Hp 

where H r is the ring Hamiltonian given in the previous section and 

n 

H p (x, y, 6i, (pi) = —k ^2 In ^1 — x sin 8j cos <j)j — y sin 9j sin <j)j — \/l — x 2 — y 2 cos Oj^j , 
j=l 

is the Hamiltonian responsible for the interaction of the pole and the ring. 
In this case 



= H - C \ J2 cosd i + ^V 1 ~x 2 -y 2 
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and the relative equilibrium at x = y = 0, 6j = 8q, 4>j = 2nj/n has angular velocity 

(n — 1) cos #0 + k(1 + cos 6>o) 
sin 2 6> 

since ^i( Xe ) = _ (n-i)co8eo+K(i+cosg )-gsin^o must vanish _ 

The second derivatives at the relative equilibrium of H can be derived from those for H r given in 
the previous section, together with: 



d 2 H p 


K 


d 2 H p 


riK d 2 H p 


aej 


1— cos 8q 


dx 2 


2 dy 2 


d 2 H p 


k cos(27rj/n) 


d 2 H p 


k sin(27r;/n) 


dxddj 


1— cos 0o 


dyddj 


1— COS 0Q 


a 2 H p 


k sin 9q sin(27rj/n) 


d 2 H p 


Ksin0o cos(27rj'/n) 


dxd(f>j 


1— COS 0Q 


dyd<i>j 


1— COS 0Q 



while the other second derivatives all vanish. Here we have used that ^ cos 2 (27rj'/n) = n/2 for 
n > 2, but for n = 2 the sum is 2. Thus for n = 2, one obtains 

d 2 g p __ 2k 9 2 Hp _ 2k cos o 

ftr 2 1— cos0o dy 2 1— cos0o' 

The following proposition gives the symmetry adapted basis for C nv (R,p) relative equilibria. 
Proposition 5.1 Let n > 3 and /i / 0. With respect to the following basis for the symplectic slice: 
(e u e 2 ,e 3 ,e i ,{a%\af,$\$ ) \2<t<[n/2}y^ 

ncos(26> )/(2fi;) 5y n 
n cos(26*o) / (2k) 8x n 

the Hessian d 2 H^\_\f(x e ) block diagonalizes into lxl blocks and two 2x2 blocks, and Lj\f block 
diagonalizes into 2x2 blocks and one 4x4 block. 

Proof. The proof is similar to the proof of Proposition 14. II □ 

These block diagonalizations enable us to prove the following stability theorem for n > 4, illus- 
trated by Figures I5^2l and l5l71 The cases n = 2 and 3 are treated afterwards. 

Theorem 5.2 A C nv (R,p) relative equilibrium with n > 4 and p ^ 

(i) is spectrally unstable if and only if 

k < kq or 8an > (n sin 2 9$ + 4(n — 1) cos 9q) 2 , 

( ii) is Lyapunov stable if 

K > Kq and clk(k + n cos Oq){k — k{) < 0, 



e\ = cos 6q Pq — sin Oq — 

e2 = sin 6*o a.g + cos 9q 

e3 = cos 9q + sin 6q f3y — 

e4 = sin 6*o — cos 9q off 
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where 



o = (n cos 6*o — n + 2)(1 + cos 6>o) 2 
ki = (n — 1) cos Oq (n sin 2 $o + 2(n — 1) cos #o)/ a 
= (cn - (n - 1)(1 + cos 2 0„))/(l + cos ) 2 



K 



J n 2 /4 i/ n is even, 

° n ~ \(n 2 -l)/4 if n is odd. 

Lyapunov stability. Following the beginning of the proof of Theorem l4.2l 
1511 that 

d 2 H(\x(x e ) = diag(A A, D) 
where D = diag({A^ xf , X^\ xf \2<l< [n/2]}*), 



Proof. We first study the Lyapunc 
we obtain from Proposition 15 . 1 1 that 



A 



911 912 

912 922 



and 



9n 

912 
922 



d 2 H s (x e ) 
d 2 H ( (x e ) 
d 2 H^(x e ) 
d 2 H/:(x e ) 
d 2 H i .{x e ) 



(a 

H »« 

(ei,ei) 
(ei,e 2 ) 
(e2,e 2 ). 



Note that D exists only for n > 4. From the previous section one has A 
additional computations give 



(0 



n£(n — £)/2 and some 



A 



2 sin 2 6» 



[-(* - l)(n - £ - 1) + (n - 1) cos 2 O + «(1 + cos O ) 2 



The eigenvalues Ar are all positive, thus D is definite if — (£ — l)(n 
k(1 + cos6» ) 2 > for all i = 2, ... , [n/2], that is if k > (([n/2] 
1) cos 2 #o)/(l + cos #o) 2 which corresponds to k > ko- 

The relative equilibrium is therefore Lyapunov stable if A is positive definite, that is if (711(722 
q 2 2 > and 911 + (722 > 0. Some lengthy computations give 



I - 1) + (n - 1) cos 2 O + 
l)(n - [n/2] - 1) - (n - 



9H922 



912 
922 



■(« - «2) 



2(l+cos 6» )- 



where a, «i are given in the theorem and k 2 = — 2(n — 1) cos 2 #o/(l + cos $o) 2 - Now we show that if 
9n 922 - 9i2 > and « > «0» then 911 + 922 > 0: we have 922(911 + 922) > since 911922 - 9i2 > °> 
and 922 > since k > ko > K2, hence 911 + 922 > 0. We proved therefore that C nv (R,p) is 
Lyapunov stable if k > kq and an(n + n cos 0q)(k — k\) < 0. 

We now study the spectral stability of the relative equilibrium. It follows from Proposition 15. H and 
the block diagonalization of d 2 H^\j^(x e ) that 




-A 







-A 



(0 



2<£< [n/2] 
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where the blocks are given up to a positive scalar factor and 



A L 



( o 





V 



-a 
-c 



) a b \ 

) c d 

6 

d / 



a 
b 
c 
d 



and 



a 

P 
7 



co (ei, es) = n cos #0 sin 2 i 
to (e2, 64) = n cos #o sin 2 i 



PQ12 

«9i2 
«922 

,o2 



7912 
7922 

7911 
7912 



o -n 2 cos 2 (26'o)/(4k) 



a; (ei, e 4 ) = w (e 2 , e 3 ) = nsm0 Q /2 
The eigenvalues (up to a positive factor) of La/ are therefore 

1 



V2 



±y/v + cr, < dbi 



2 < £ < [n/2] 



where ^ = a 4 + 4a 2 6c — 2a 2 d 2 + 46ai 2 + <i 4 + 8acZ6c and cr 



-a 2 — 26c — d 2 . The eigenvalues 

Xg'X£ J are all purely imaginary if and only if k > kq. After some lengthy but straightforward 
computations we obtain that 

v = ^(l +u ) 2 (2n 2 -l) 8 ( K + nn) 4 (l-n) 2 



x 



[—8(1 + u) 2 (nu + 2 - ti)k + (nu 2 - 4(n - l)u - n) 2 



a 



-Jfe(2u 2 -1)*(k + nu) 2 

x [—4(1 + u) 2 (nu + 2 — n)K + n 2 u 4 — 4n(n — l)u 3 + 
+2(3n 2 - 8n + 4)n 2 + 4n(n - l)u + n 2 ] 

where u = cos 9q. One can check that if v > 0, then yfu + a < and the eigenvalues are purely imag- 
inary. If v < 0, then the eigenvalues have a non-zero real part. Thus the eigenvalues ±y±^/z/ + a 
are purely imaginary if and only if v > which is equivalent to 8cjk < (n sin 2 6*o + 4(n — 1) cos #o) 2 - 

□ 



A spectrally stable relative equilibrium for which the Hessian d 2 H^\j^f(x e ) is not definite is said to 
be elliptic. Note that in principle an elliptic relative equilibrium may be Lyapunov stable, but if there 
are more than 4 vortices then it is expected to be unstable as a result of Arnold diffusion. Moreover 
an elliptic relative equilibrium typically becomes linearly unstable when some dissipation is added to 
the system [DR021. 

Corollary 5.3 A C nv (R,p) relative equilibrium with n > 4 and [i ^ is elliptic if and only if 
k > kq, ok(k + ncos 0q)(k — K\) > and 8an < (nsin 2 #o + 4(n — 1) cos #o) 2 - 

Discussion of results for n > 4 See Figure l5~T1 

• If the sign of the vorticity of the polar vortex is opposite to that of the ring then there are stable 
configurations with k < only for n < 6. Conversely configurations with n < 6 and 9q close 
to 7T, ie with the ring close to the opposite pole, are Lyapunov stable for all k < 0. 
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Figure 5.1: Bifurcation diagrams for C nv (R,p) relative equilibria. The bifurcation diagrams for 
n > 7 are similar to that for n = 8, while those for n = 4, 6 are similar to that for n = 5. The circles 
represent the eigenvalues of the mode 1, while the cross represent those of the mode [n/2]. A relative 
equilibrium can be Lyapunov stable (S), linearly unstable (U) or elliptic (E). 
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• The region of Lyapunov stability is larger when the vorticities of the pole and the ring have the 
same sign (k > 0). The stability frontiers in the upper-left corners of Figure 15. II go the infinity 
when #o g° es to arccos(l — 2/n). It follows that for n > 4 and 9q > arccos(l — 2/n), the 
relative equilibria are Lyapunov stable for all sufficiently large n. Thus, a ring of vortices is 
stabilized by a polar vortex with a sufficiently large vorticity of the same sign as the vortices in 
the ring. Note that for 4 < n < 6 and k positive, but sufficiently small, a ring near the opposite 
pole is only elliptic and may not be Lyapunov stable. 

• The limiting stability results for 9q = 0, ie when the ring is close to the polar vortex, coincide 
with the stability of a planar n-ring plus a central vortex, see [CS99 ] and [LP|. This is also true 
for n = 2 and n = 3. 

• One of the main stability boundaries corresponds to the mode [n / 2] and is analogous to the sta- 
bility boundary for a single ring. When n is even stability is probably lost through a pitchfork 
bifurcation to a relative equilibrium of type Cil v (R, R',p) consisting of two staggered |-rings 
and a pole as (k, 9q) passes through this boundary. When n is odd there is an analogous tran- 
scritical bifurcation to relative equilibria with only a single reflectional symmetry which fixes 
two vortices and permutes the others. These are denoted by C/ l ( 2 ^i?, 2E) in [LMR01 1. A 
nice illustration in the case n = 3 can be found in Figure 8 of [CMS03 1. 

• The other stability boundary corresponds to the mode 1. Stability is lost through a Hamiltonian- 
Hopf bifurcation: two pairs of imaginary eigenvalues 'collide' and leave the imaginary axis. For 
a detailed description of the bifurcations that can be expected in this case see [ vM85 1. 

• Note that it also happens that pairs of £ = 1 eigenvalues pass through zero without leaving the 
imaginary axis. In this case the relative equilibrium changes from being Lyapunov stable to 
elliptic or vice versa and these stability changes are accompanied by bifurcations. 

• Finally we note that when k crosses zero eigenvalues change sign without crossing zero due to 
the fact that the symplectic form becomes degenerate for k = 0. 

Discussion of the case n = 3 See Figure l5~2l 

We assume in this discussion that the relative equilibria have non-zero momentum. For n = 3, by 
the proof of Proposition 15 .21 we have d 2 H^\j^(x e ) = diag(yl, A) and Ljy = A^. Hence C3 V (R,p) is 
Lyapunov stable if ok(k + 3 cos 9q)(k — k\) < 0, and spectrally unstable if and only if 

8an > (3 sin 2 9 + 8 cos 9 ) 2 , 

where a = (3 cos 9q — 1)(1 + cos 9q) 2 as in Theorem l5.2l These results are illustrated in Figure |5"21 
Notice that a polar vortex destabilizes a 3-ring if either the polar vortex is in the same hemisphere as 
the ring and has a sufficiently strong vorticity of the same sign as the ring, or the polar vortex has 
the opposite sign vorticity and the ring lies in an interval containing 9q = 2n/3 that grows as the 
magnitude of the polar vorticity increases. Outside these regions there is a patchwork of regimes in 
which the relative equilibrium is either Lyapunov stable or elliptic. 

Discussion of the case n = 2 See Figure |5~21 

The C2v(R,p) relative equilibria are isosceles triangles lying on a great circle, and for 9o = 2ir/3 the 
triangle becomes equilateral. We again discuss the stability of those with non-zero momenta. Indeed, 




Figure 5.2: Bifurcation diagrams for C2 V (R,p) and C^ v (R,p) relative equilibria. The circles repre- 
sent the eigenvalues of the mode 1 (some are filled in order to trace them). A relative equilibrium can 
be Lyapunov stable (S), linearly unstable (U) or elliptic (E). 
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any 3-vortex configuration with zero momentum is a Lyapunov stable (relative to SO(3)) relative 
equilibrium since the reduced space is just a point IPa92l . For n = 2 the symmetry adapted basis is 
(k 56\ — k 502 — 2 cos #o 5x, k 5(j>\ — k d(f>2 — 2 sin 9q Sy). Following the proof of Theorem 15 .21 we 
obtain after some straightforward computations that C2 V (R,p) is Lyapunov stable if 

(1 + 2cos6»o)[(1 + cos6>o) 2 k + cos<9 (2 + 3cos(9o)] < 0, 

and linearly unstable if the inequality is reversed. See Figure l5T2l 

• There are two stable regions. For 9q < 27r/3 the relative equilibria are stable provided the polar 
vorticity is less than a certain 9 dependent critical value, while for 9$ > 2ir/3 they are stable 
for all polar vorticities greater than a critical value. As 9q — > n this value goes to — oo. 

• For 6*o = 7r/2, where the 2-ring is equatorial and the isosceles triangle is right-angled, they are 
stable if and only if k < 0. This is in agreement with [PM98 Theorem III.3], with Ti = T2 = 1, 
and T3 = k. 

• The restricted three vortex problem The range of stability when k = does not coincide with 
the range of stability for a single ring. Indeed the C2 V (R) relative equilibria are Lyapunov stable 
for all co-latitudes (see Theorem l4.2t while C2 V (R,p) is unstable for k = and 9q S (0, 7r/2). 
This means that if we place a passive tracer (or ghost vortex) at the North pole and a ring of 
two vortices in the Northern hemisphere, then the passive tracer will be unstable. 

Remark 5.4 The stability of C nv (R,p) relative equilibria has also been studied in [CMS03|. How- 
ever our method differs significantly from theirs in that we consider the definiteness of the Hessian 
d 2 H^\j^(x e ) on the 2n dimensional symplectic slice, while in [CMS03] the authors determine condi- 
tions for the Hessian to be definite on the whole 2(n + 1) dimensional tangent space. The result is 
that we prove the relative equilibria to be Lyapunov stable in a larger region of the parameter space. 
Notice in particular that for n < 6 our results say that a positive vorticity n-ring near the south pole 
is Lyapunov stable if the north pole has either negative or sufficiently positive vorticity. However in 
[ CMS03 1 only the case of negative north polar vorticity is shown to be Lyapunov stable. In this paper 
we also give criteria for when the relative equilibria are unstable by considering the eigenvalues of the 
linearisation L/v. 

6 Stability of a ring and two polar vortices: C nv (R, 2p) 

In this section we consider a relative equilibrium x e of symmetry type C nv (R, 2p), that is configura- 
tions formed of a ring of n vortices of strength 1, together with two polar vortices p^, ps of strengths 
kn, ks respectively at the North and South poles. We assume without loss of generality that the ring 
lies in the Northern hemisphere. 

We obtain analytic (in)stability criteria for the relative equilibria with respect to the £ > 2 modes, 
which of course give sufficient conditions for genuine instability. A numerical investigation is needed 
for the £ = 1 mode and hence to provide stability criteria; this is being pursued separately. 

The Hamiltonian is given by 

H = H r + H PN + H ps + Hns 
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where H r is given in Section|4]and 

H PN = ~K N YZ=1 Ki ln ( 1 ~~ sin ®i COS & XN ~ sin 0i Sin ~ vA ~~ X N ~ Vis 



COS ( 



H ps = -K S ^ =1 Kiln(l - smOiCosfa x s - smOisinfa y s + yl - x 2 s - ygCosOi) 
H N s = -k n k s ln(l - x N x s - y N y s + yjl - x 2 N - y 2 N ^l - x 2 s - y 2 s ), 

and the augmented Hamiltonian is: 



= H - £ ^2 ^ cos 6i + K N y 1 - x 2 N - y 2 N - k S \J 1 - x 2 s - y 2 s 



\i=l 



The angular velocity of the relative equilibrium at xn = UN = %s = Us = 0, 8j = 6o, <ftj = 
2-irj/n has angular velocity 

(re — 1) cos #o + «jv(1 + cos #o) — K s(l — cos 0o ) 
sin 2 #o 

The second derivatives of H at the relative equilibrium can be derived from those for H r (Section 
131), those for H p (Section|5ll, together with: 

d 2 H NS d 2 H NS d 2 H NS d 2 H NS d 2 H NS d 2 H NS 



dx 2 N dy 2 N dx 2 s dy 2 s dx N dx s dy N dy s KnKs ^ 2 

while the other second derivatives of Hns vanish. 

As in the previous sections, we can choose a symmetry adapted basis of the symplectic slice such 
that the matrices d 2 H^\j\f(x e ) and L/v block diagonalize. 

Proposition 6.1 Let n > 3 and /J, ^ 0. In the following basis for the symplectic slice, 

(e 1 ,e 2 ,e 3 ,e 4 ,e 5 ,e 6 ,{af\a i l\$\$ ) I 2 < £ < [re/2]}*) 

where 

y - ncos(26 )/(2K N ) 5y N 

^ -ncos(20 o )/(2K S )Sy s 
(i) 

y - ncos(2d )/(2K N ) 5x N 

^ -ncos(20 o )/(2K S )Sx s 
(i) 

the Hessian d 2 H^\j^(x e ) block diagonalizes into lxl blocks and two 3x3 blocks, and Lj\f block 
diagonalizes into 2x2 blocks and one 6x6 block. 

Proof. The proof is similar to that for a single ring (see Section EOl and Proposition l4.lt . □ 

The mode 1 = 1 gives a 3 x 3 block from which, unfortunately, we can not derive a useful 
formula for stability analogous to that for a single polar vortex. However, we can derive formulae for 
the stability of the other modes, and thereby obtain the following sufficient condition for instability, 
illustrated by Figures l6~Tl and l6~2l 
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Theorem 6.2 A C nv (R, 2p) relative equilibrium with n > 4 and u ^ is linearly unstable if 
kjv(1 + cos #o) 2 + — cos 6q) 2 < c n — (n — 1)(1 + cos 2 #o)> 

where 

( n 2 /4 if n is even, 

Cn ~\(n 2 -l)/4 if n is odd, 

and is stable with respect to the £ > 2 modes if this inequality is reversed. 



-co 

80- 




Figure 6.1: The relative equilibria C nv (R, 2p) in the Northern hemisphere are unstable 'below' this 
ruled surface in (6q, ks, KAr)-space, shown in the figure for n = 4. Above the surface the relative 
equilibrium is stable with respect to all the £ > 2 modes. 



Proof. The proof is similar to that for Theorem l5.2l Following the notations of the proof of Theorem 
K2l we have A f = ni{n - I) /2 and 

\i e) = \-(£ - l)(n - £ - I) + (n - 1) cos 2 6 + k n (1 + cosOo) 2 + k s (1 - cosOo) 2 } . 

2 sin 9q 

The relative equilibrium is linearly unstable if there exists £ > 2, such that \\ < 0. Since the highest 
Xy is for £ = [n/2], the relative equilibrium is linearly unstable if 

-fln/2] - l)(n- [n/2] - 1) + (n - 1) cos 2 6> + k n (1 + cos 6» ) 2 + k s (1 - cos 6» ) 2 < 0, 

and is stable with respect to the £ > 2 modes if this inequality is reversed. This gives the desired 
criterion. □ 
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Stability of the I > 2 modes From the theorem we can deduce the following results about the 
(in)stability of the C nv (R, 2p) relative equilibria with respect to the I > 2 modes. These modes only 
occur for n > 4. We continue to assume the ring lies in the Northern hemisphere. 

• In the limiting case as the ring converges to the North pole (9q = 0), for all values of k$ the 
relative equilibria are linearly unstable if kn < \{c n — 2n + 2). This agrees with the instability 
of a ring and single pole when 'k < in Proposition 15 .21 

• At the opposite extreme, when the ring is at the equator (6q = tt/2) they are linearly unstable 
if ktv + ks < c n — n + 1. The right hand side of this inequality is non-negative for all positive 
integers n, and so the 'equatorial' C nv (R, 2p) relative equilibria are unstable if the total polar 
vorticity has opposite sign to that of the ring. If kn + k-s > then the critical ratio of the total 
polar vorticity to the total ring vorticity needed to stabilize the t > 2 modes grows linearly with 
n. 

• For all n > 4 the relative equilibria are unstable for all latitudes in the Northern hemisphere if 
kn < \ (c n — 2n + 2) and kn + k$ < c n — n + 1. In particular, for n > 7 the relative equilibria 
are unstable for all 6q if kn < and k$ < 0. 




Figure 6.2: (a) The shaded regions depict the values of the polar vorticities for which all the relative 
equilibria in the Northern hemisphere are unstable: the dark region is for n = 4, the next for n = 8 
and the lightest for n = 10. (b) Demonstrates that above each shaded region of (a) the corresponding 
relative equilibria near the North pole are unstable, while to the right it is the relative equilibria near 
the equator which are unstable. The instabilities are all due to the I > 2 modes. 

To determine whether there are in fact stable relative equilibria it is necessary to evaluate the 
eigenvalues arising from the £ = 1 mode. This is work in progress, and preliminary numerical 
investigation suggests: 

• For all n and for all sufficiently large and positive polar vorticities there are ranges of 9q with 
elliptic relative equilibria; 

• For all n and for kn sufficiently positive and ks < there are Lyapounov stable relative 
equilibria in the Northern hemisphere. 
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7 Stability of two aligned rings: C nv (2R) 

In this section we consider relative equilibria x e of symmetry type C nv (2R), that is configurations 
formed of two 'aligned' rings of n vortices each. We can assume without loss of generality that the 
vorticities of the vortices in the first and second ring are 1 and k, respectively, and we denote their 
co-latitudes by 9q and 9\ . We can also assume that the ring of vorticity 1 and co-latitude 9q lies in the 
Northern hemisphere, 6q £ (0, n/2]. The first question to answer is for which values of the parameters 
(9q,9i,k) is the configuration C nv (2R) a relative equilibrium? It was shown in [LMR01| (p. 126) 
that for given k > and each fi with \u,\ < n\l + k\ there is at least one solution for (9q, 0{) with 
n cos #o + n-K cos 9± = p and with 9q < #i an d at least one with 9\ < 9q. We now make this more 
precise. 

The fixed point set Fix(Ga; e ) is parameuized by x := cos^o and y := cos0\. Denote by F the 
restriction of a function F to Fix(G Xe ). It is straightforward to show that the Hamiltonian can be split 
in such a way that 

H = H n + kH 12 + n 2 H 22 

where Hn, H\ 2 , H 22 do not depend on k, H%i does not depend on y and H 22 does not depend on x. 
The following proposition shows that for almost every pair (9q, 9\) there exists a unique n such that 
the C nv (2R) configuration with parameters (6>o, 9\, k) is a relative equilibrium. 

Proposition 7.1 Let x e be a C nv (2R) configuration with parameters (9q, 9%, k). 

1. There exists a unique K G M* such that x e is a relative equilibrium if and only if both the 
following conditions hold: 

(dH 12 dH n \ ft ,, n (9H 22 dH 12 \ 

\W - "to" ) ( COS ^' C ° S ^) ^ °' [-Qy- ~ "ftT ) (^^0,COS^) * 0. 

2. The configuration x e is a relative equilibrium for all k£1* in the degenerate case when both 
the following conditions hold: 

(dH l2 dH n \ (dH 22 8H 12 \ 

I —q^ q^- I (cos 0q, cos 9\ ) = 0, I — — I (cos 9q , cos 9\) = 0. 

3. In both cases the angular velocity £ of x e satisfies 



£ = 1 dHll (x e ) I r dHl2 (x 
n \ dx 6 dx 



In particular, if in addition d x Hn(x e ) and d x H\ 2 (x e ) are non-zero, then there exists a unique 
k£M* such that £ = and so x e is an equilibrium. 

Proof. Since H + is a G Xe -invariant function (see Section |2j) the Principle of Symmetric Crit- 
icality [P79| implies that x e is a relative equilibrium if and only if it is a critical point of H + £5>. It 
follows from $ = n(x + ny) that d(H + £$)(a; e ) = is equivalent to the pair of equations: 

(dH 22 dH 12 .\,dHu dH lx 

~dy- {Xe) ~ ^X~ {X&) + ^y~ {Xe) ~ ^x~ {Xe) = ° 
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1 (dHn dH 12 \ 

5 = - (^e) + K— — (x e ) . 

n \ ax ox / 

The proposition follows easily from these. □ 

For example, in the case n = 4 the degenerate case occurs when the two rings form the vertices of 
a cube. Hence for any values of the vorticities of the two rings the "cube configuration" is a relative 
equilibrium. However among this family of relative equilibria only one is an equilibrium, namely 
the one for which the two rings have the same vorticities, k = 1, which corresponds to the Oh{f) 
equilibrium |LMR01 1, a cube formed of identical vortices. 

For 6*i = 7T — 6*o, the configuration has an extra symmetry and its symmetry type is D n f l (2R). 
Such a configuration is a relative equilibrium if k = — 1, the two rings have opposite vorticities. The 
existence and stability of such relative equilibria were studied in ILP02I . 

With the help of the discussion of Section 1331 we performed a G Xe -invariant isotypic decom- 
position and found that the symmetry adapted basis for the symplectic slice at a C nv (2R) relative 
equilibrium with n > 3 and fx ^ is 

(ei, e 2 , e 3 , e 4 , e 5 , e 6 , e 7 , e 8 , {B e | 2 < £ < (n - l)/2}) 

for n odd, while it is 

( (T3 I 1 ^ o ^ /o 11 (n/2) (n/2) (n/2) (n/2)\ 

I ei,e 2 ,e 3 ,e4,e5,e6, ,e 7 ,e 8 ,{^ | 2 < £ < n/2 - 1} ',a^ , ) 

for n even, where 

_ (o) _ (o) 

ei ~~ a O,0 a l,0 

• a (0) ■ /i (0) 
e 2 = ac sin t/i - sm do a\ e 

e 3 = sin sin 9\ (k cos 9\ ol^'q + cos 8q a^g) 

+ cos O cos 9 1 (k sin 0i 0^ + sin O 
e 4 = k cos 0i a Q q — cos 0o e 
e 5 = k sin 0i - sin O 0^ 
e§ = sin 0o sin0i (kcos0i 0^1 +cos0o 0ig) 

— cos 0o cos 0i (k sin 0i a^p, + sin 0o ct^) 

e-J = K COS 01 0Qg — COS 00 0^ g 

■ad) -a (i) 
e 8 = Ksmffx a> Q ^ - sm0 o 

and 

r> _ f W _ W fl W _j_ M , (0 fl W _ ail) 

&l — | a 5 e 0*1,0) Po,0 "r ^1,0' a O,0 "r a \,4>i l J 0,6 P\fii 

(0 , CO _ M _ (-0 , \ 
a o,e + a i,6»' Po,(/> Pi,</>i a o,</> a i,0> ^o,e + "1,0 j • 

The adapted basis for n = 2 is simply (ei, e 2 , e 4 , e 8 ). 

Remark. Almost all C n „(2i?) relative equilibria have a non-zero momentum. Indeed \i = iff 
x + m/ = 0, and from the expression of k one can show that this last equation is the equation defining 
an algebraic curve in variables (x,y) G [0, 1) x (—1, 1) ~ FixC n „. 
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Figure 7. 1: Sign of k and stability frontier for C nv (2R) relative equilibria. The degenerate case occurs 
when the curves k = and k^ 1 = intersect each other. The figure is for n = 6, but is similar for 
other values of n (concerning the sign of n). 



With respect to this basis d 2 H^\j^(x e ) block diagonalises into two 4x4 blocks for each I = 
1 . . . [(n — l)/2]), together with two 2x2 blocks for £ = n/2 when n is even. The linearisation La/ 
block diagonalises into half as many blocks of twice the size. In order to calculate the stability of the 
relative equilibria, we ran a Maple programme to compute numerically the eigenvalues of each of 
the blocks of d 2 H^\j\f(x e ) and Lj^. The results are summarized for n = 2 ... 6 in Figure l7~2l Figure 
l7.1l shows how the sign of k varies for relative equilibria with different values of 9q and 9\. 

Discussion of results 

• The numerical results suggest strongly that the relative equilibria C nv (2R) are never stable if 
the two rings lie in the same hemisphere (Figure 17.21) or have the same sign vorticity (Figure 



• The stable configurations are for 9$ close to and Q\ close to tt, so the ring of vorticity 1 is 
close to the North pole and the other ring is close to the South pole with a vorticity close to —1. 
Thus the two rings 'look like' two polar vortices. It is well known that any configuration of two 
vortices is Lyapunov stable. 

• As n increases the region of stability decreases in size. Numerical experiments with n > 7 
suggest that in these cases the relative equilibria are never stable. 

• For n = 2, 4 and 6 stability is first lost by a pair of imaginary eigenvalues of the £ = n/2 
block of Lj\f passing through and becoming real. For n = 3 and 5 close to the D n h(2R) 
relative equilibria a pair of imaginary eigenvalues of the I = (n — l)/2 block passes through 
but remains on the imaginary axis, so the stability changes from Lyapunov to elliptic. This 
imaginary pair then collides with another pair, and all move off the imaginary axis to form a 



mi . 
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Figure 7.2: Stability results for C nv (2R) relative equilibria. The curve plotted for each n is the 
stability frontier: on one side the relative equilibria are Lyapunov stable (S), while on the other side 
they are linearly unstable (U), with a gap of elliptic stability (E) between the two in the odd case. For 
n > 7, it seems likely that the relative equilibria are all unstable. 



complex quadruple and create instability. It seems likely that this behaviour also occurs away 
from the D n h(2R) relative equilibria, but in a region too small to be seen in the figure. 



8 Stability of two staggered rings: C nv (R,R') 

In this section we consider relative equilibria formed of two rings of n vortices each of strengths 1 
and k and co-latitude 9q and 6\ respectively. They differ from those of the previous section in that 
the rings here are "staggered", that is they rotated relative to each other with an offset of n/n. Their 
symmetry type is C nv (R, R'). As in the previous section we can assume without loss of generality 
that the ring of vorticity 1 and co-latitude 9q nes m the Northern hemisphere. 

An analogue of Proposition 17. 1 1 also holds for C nv (R, R') configurations: for almost every pair 
(#0)#i) there exists a unique k such that the corresponding C nv (R, R') configuration is a relative 
equilibrium. With the notation of the previous section, in the non-degenerate case the angular velocity 
and k satisfy 

k = -(d x Hn - dyH 12 )J(dyH22 - d x H 12 )(x e ), 
£ = n - l (d x Hu + Kd x H 12 )(x e ). 
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There exist also degenerate cases. When 9\ = 9q the configuration forms a single ring with 2n vortices 
with k = 1: all the vortices have the same vorticity. These are the relative equilibria of type C2nv(R) 
studied in Section l3~2l For 9\ = ir — 9q, the configuration has an extra symmetry and its symmetry 
type is D n d(R, R'). In this case k = — 1, the two rings have opposite vorticities. The existence and 
stability of such relative equilibria were studied in [LP02|. 

With the help of the discussion of Section 13.31 we found that a symmetry adapted basis for the 
symplectic slice at a C nv (R, R') relative equilibrium with n > 3 and fj, ^ is given by: 



where the expressions of e\,...,e% and Bi remain as in the previous section. The corresponding 
symmetry adapted basis for n = 2 is simply (ei, e§). As in the previous section, we can show 

that almost all C nv (R, R') relative equilibria have non-zero momenta. 

As for the aligned rings, we ran a MAPLE programme to determine the stability of the relative 
equilibria. The results are summarized in Figure UTTl for n = 2 . . . 6. 

Discussion of results 

• Numerical experiments suggest that stable relative equilibria only exist for n < 6. 

• For n = 5 and 6 the relative equilibria C nv (R, R') are stable only if the two rings lie in the 
same hemisphere but are sufficiently far apart. 

• For n < 4 these stable regions extend to include relative equilibria with the rings in different 
hemispheres. However, contrary to the case C nv (2R^, the stable regions are far from the line 
9\ = 7r — 9q corresponding to D n d(R, R') relative equilibria. 

• For n = 2 and 3 there is also a stable region with the two rings in the same hemisphere and 
close to each other. This includes the stable C^i?) and Cq v (R) relative equilibria discussed 
in Section l3~2l 

• Note also that for n < 6, there exist stable relative equilibria (for some values of k) in any 
neighbourhood of (#o, #i) = (0, 0), that is with the two rings close to the North pole. 

• A study of the sign of k shows that when 3 < n < 6 the relative equilibria with k < are 
all unstable. However for n = 2 there exist relative equilibria with k positive in the stable 
region corresponding to the two rings both being relatively close to the equator, but in opposite 
hemispheres. 

• The tetrahedral equilibrium with all 4 vortices identical is Lyapounov stable, and in Fig. 18.11 
(with n = 2) lies at the point where the two stable regions meet the two unstable regions on 
the D 2 d(Ri i?')-locus. The analogous point with n = 3 corresponds to the stable equilibrium 
consisting of 6 identical vortices lying at the vertices of an octahedron Oh(v). 



(ei,e2,e 3 ,e4,e5,e 6 ,e 7 ,e 8 ,{^ | 2 < t < (n - l)/2}) 



for n odd, while for n even one is given by: 



(ei,e2,e 3 ,e4,e5,e 6 ,e 7 ,e 8 ,{^ | 2 < t < n/2 - 1} 
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Figure 8.1: Stability results for C nv (R, R') relative equilibria. The curves plotted are stability fron- 
tiers: on one side the relative equilibria are Lyapunov stable (S), while on the other side they are 
linearly unstable (U). For n > 7, it seems likely that the relative equilibria are all unstable. 
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